Anomalous thermodynamics at the micro-scale 
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Particle motion at the micro-scale is an incessant tug-of-war between thermal fluctuations and 
applied forces on one side, and the strong resistance exerted by fluid viscosity on the other. Friction 
is so strong that completely neglecting inertia - the overdamped approximation - gives an excellent 
effective description of the actual particle mechanics. In sharp contrast with this result, here we 
show that the overdamped approximation dramatically fails when thermodynamic quantities such 
as the entropy production in the environment are considered, in presence of temperature gradients. 
In the limit of vanishingly small, yet finite inertia, we find that the entropy production is dominated 
by a contribution that is anomalous, i.e. has no counterpart in the overdamped approximation. 
This phenomenon, that we call entropic anomaly, is due to a symmetry-breaking that occurs when 
moving to the small, finite inertia limit. Anomalous entropy production is traced back to futile 
phase-space cyclic trajectories displaying a fast downgradient sweep followed by a slow upgradient 
return to the original position. 



Life at the micro-scale flows under one law: the fluid 
gives and the fluid takes away. This arbitrary tyrant 
lavishly bestows energy to suspended particles through 
molecular collisions while incessantly draining it from 
them through friction. The result is the erratic mov- 
ment of microscopic particles that goes under the name 
of Brownian motion. 

The theory of Brownian motion was developed by Ein- 
stein, Smoluchowski and Langevin a little over a century 
ago pQ . A central result of this theory is the overdamped 
approximation, which says that inertia can be ignored if 
mass is small, or friction is large. The motion of a Brow- 
nian particle, which obeys Newtonian mechanics and is 
driven by collisions and external forces, is thus reduced 
to a first-order diffusion equation. The overdamped ap- 
proximation successfully describes the mechanics on the 
microscale and is very widely used [2j |3] . 

With the advent of micromanipulation it has become 
possible to measure and control the positions of indi- 
vidual Brownian particles and other small systems [4]- 
[7]. Thermodynamic concepts such as heat, work and en- 
tropy production have hence taken a meaning for single 
systems [8, 9 . These developments are the foundation 
of the new emerging field of stochastic thermodynamics, 
where the most striking results obtained to date are fluc- 
tuation relations, recently reviewed in, e.g., Refs. [T0UT2] . 
Here we show that in this setting the overdamped approx- 
imation fails, as soon as the temperature field varies in 
space. 

Indeed, although the overdamped approximation cor- 
rectly yields the trajectories of the Brownian particles 



in space, it incorrectly estimates the entropy production. 
This failure is traced back to seemingly innocent correc- 
tions that, while having a negligible impact on trajecto- 
ries, eventually dominate the entropy production in the 
long run, even in the limit when inertial effects go to zero. 
As a result, fluctuation relations themselves take a very 
nontrivial form in the limit of small yet finite inertia. 

We dub this phenomenon entropic anomaly, in anal- 
ogy with similar anomalies encountered in physics. The 
best known example in classical physics is the viscous 
dissipative anomaly. The energy dissipation in a fluid 
flow remains finite even in the limit of arbitrarily small 
viscosity, whereas it vanishes when viscosity is exactly 
zero [13 . The viscous anomaly reflects the loss of time- 
reversal symmetry of fluid dynamics when going from 
the inviscid case to the viscous one. Quantum anoma- 
lies arise when a current conserved at the classical level 
(Planck's constant set to zero) is not conserved anymore 
at the quantum level (see, e.g., Ref. [14] for a general dis- 
cussion and Ref. [15 for a simple example from molec- 
ular physics). We show here that the entropic anomaly 
is associated to the breaking of a symmetry of the zero- 
inertia, overdamped dynamics, when small finite inertia 
is considered. The symmetry is the joint reversal of time 
and particle velocity. 

To illustrate our result, let us focus on a minimal 
thought experiment. A microscopic particle is suspended 
in a fluid inside a vessel with reflecting walls. The two 
opposite sides of the vessel are in contact with heat reser- 
voirs at different temperatures. The fluid is motionless 
and displays by a static, smooth temperature profile T(x) 
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(e.g. a linear one). For the sake of simplicity, no exter- 
nal force is applied to the particle and the friction co- 
efficient is taken independent of particle position. We 
emphasize that the assumptions on temperature profiles, 
friction coefficients and absence of external forces are not 
restrictive, and refer the discussion below and the Sup- 
plementary Material for a general discussion. 

The motion of the suspended particle is then governed 
by the Langevin-Kramers equations for position X t and 
velocity V t [16] 

x t = v t 

V t = - 7 Vt + y/2T(X t hTh 1 ' 

where r} t is a Gaussian, zero-mean white noise, i.e. 

(44) = s»s(t-t'). 

The heat released by the particle to the fluid along a 
trajectory from time t' to time t reads Q = — \V t \ 2 /2 + 
\V t f\ 2 /2. In the case at hand, since there is no force, no 
work is done on the particle and no potential energy is 
stored [T7] , The entropy of the particle is defined as the 
state function 

S p (x, v, t) = — \np(x, v, t) (2) 

where p is the probability density of the particle position 
and velocity at time t [8 j. In other words, p is the solu- 
tion of the Fokker-Planck equation associated to Eq. 0. 
The entropy produced by the particle in the environment 
(the fluid) along a path is the integral of the released heat 
divided by temperature 

Senv = - J^T^VrOVrdT (3) 

where o denotes the Stratonovich product (the integral 
is regularized according to the midpoint rule). The total 
entropy produced along the process is therefore 

Stot = S p (Xt, t) — S p (Xf , t') + S env (4) 

and is positive on average as prescribed by the second 
law of thermodynamics [I8j [19] . 

We now turn our attention to the overdamped dynam- 
ics. It is possible to prove by means of asymptotic tech- 
niques (see [9j [20j [21] and the Supplemental Material) 
that the the spatial trajectory X t in tends - in the 
limit of small inertia and in the probabilistic sense - to 
the solution of 

X t = y/2T{X t )h m (5) 

The correct interpretation of the equation above is that 
the product on the right-hand-side has to be taken with 
the Ito, nonanticipative convention. The velocity fol- 
lows the local Maxwell-Boltzmann distribution w(v\x) = 
(27rT(x))- 3 / 2 eM-\v\ 2 /(2T(x))}. 



The Fokker-Planck equation associated to Eq. can 
be interpreted as the mass-conservation equation for a di- 
lute colloidal suspension of noninteracting particles. The 
flux of particles is J = — T7 _1 Vp— 7 _1 pVT and features 
the contributions of the osmotic force —TV In p, and of 
the thermophoretic force — VT. For times larger than 
T 2 7/T, where L is the size of the vessel that encloses the 
fluid, the probability density of particle position reaches 
the equilibrium p eq oc T _1 with a zero-flux balance of 
osmotic and thermophoretic forces. 

Stochastic thermodynamics can be formulated for the 
overdamped dynamics 0. The particle entropy is 

S^ = -In p(x,t) (6) 

where p solves the Fokker-Planck equation associated to 
Eq. ([5]). When the particle moves from a region where 
the probability density p is high to one where it is low, 
its entropy increases. The entropy produced in the envi- 
ronment is the work done by the thermophoretic force on 
the particle along a trajectory, divided by temperature 

S^ over) = - C VT(Xr) o X dr - In ( T ( Xt 'A (7) 
env ~ J t > T(X T ) 7 \T(X t ) J [() 

and it is positive when the particle moves from a hot 
region into a cold one. When the system reaches equilib- 
rium the total entropy production 

s { tot er) = s ( p over) (x t , t) - (x t - ,t') + siz er) (8) 

tends to zero along any trajectory since these two contri- 
butions compensate exactly. The average rate of entropy 
production is 

where the integral is over the volume of the vessel. It 
vanishes in equilibrium since pT is then a constant. This 
result coincides with the expression obtained by macro- 
scopic non-equilibrium thermodynamics for the entropy 
production by a dilute particle suspension [22] . 

We now move to the description of the results of the 
present work, where the limit of small yet finite inertia 
is considered. The derivations follow standard asymp- 
totic expansion methods and are detailed in the Sup- 
plementary Material. The only assumption is that the 
length-scale of variation of the temperature be larger 
than the typical distance travelled during the friction 
time 7 _1 , i.e. the gradients must not be exceedingly 
large | VT| <C Tj/v - 7T 1 / 2 [2J3]. 

The average rate of entropy production of the 
Langevin-Kramers process 0, in the limit of vanishing 
inertia, is 
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FIG. 1. Anomalous entropy is produced during futile cy- 
cles of fast downgradient sweeping and slow upgradient climb- 
ing. Along a microscopic trajectory that starts and ends at 
the same temperature (full red line), the overdamped con- 
tribution to entropy vanishes according to eq. Q. However, 
anomalous entropy is produced because of heat absorption at 
higher temperatures, followed by heat release at smaller ones. 
The reversed cycle (blue dashed line) with negative anomalous 
entropy prod uction is much less likely to occur in agreement 
with eq. |l4). 



This exact asymptotic expression has to be contrasted 
with the overdamped approximation ([9|. It differs by 
an additional, positive contribution, which actually con- 
trols the asymptotic rate of entropy production since 
d{S\°t er ^) /dt vanishes as equilibrium is approached. This 
is the most conspicuous effect of the entropic anomaly. 

It is indeed possible to isolate the source of the anomaly 
from the total entropy production Q, S to t — S reg + 
Sanom, where 



S„ 



(5T(X r ) 



2T(X T ) 2 



V T -VT(X T )dr (11) 



and to show that it arises from the S env contribution. 
The remainder S reg has a regular limit and tends to the 
overdamped entropy (JsJ) . In the limit of small inertia the 
anomalous entropy obeys the identity 



(exp(-S anom )) = 1 
and its average rate of production is 



6 J H 7 T 



(12) 



(13) 



that gives the rightmost term in (10). It is worth point 



ing out that the anomalous contribution cannot be elim- 
inated by a suitable redefinition of the overdamped en- 
tropy. Indeed, it does not exist any sequential func- 
tional of the trajectories Q that gives the correct limit- 
ing statistics of S anom (see Supplemental Material). 

What is the interpretation of the anomaly ? Entropy 
is a measure of irreversibility. It is possible to show in 



very general terms that the entropy produced in the en- 
vironment along a trajectory quantifies the probability 
of the reversed path relative to the forward one (see [18] 
for a precise statement for diffusion processes). In the 
Langevin-Kramers case, the reversal consists in the inver- 
sion of the arrow of time and of the direction of velocity 
(see Figure 1) and one has 



P(X [t , t , ] ,-Vj t , t , ] ) 



exp ( S env ) 



P{X [t , A ,V [t , A ) 
Similarly, for the overdamped dynamics 



p(o*er) (X[t>t/]) 
p{o»er)(X [t , A ) 



exp 



/ o(over) \ 

I ^ env J 



From the ratio of these two identities it follows that, 
in the limit of vanishing inertia, the anomalous entropy 
gives the relative weight of conditional probabilities 



P(Vj t ,, f] |X [f ,, f] ) 



exp ( Sanom) 



(14) 



This result shows that the anomaly emerges from the 
breaking of time and velocity reversal symmetry at a 
given spatial position. 

The anomaly can be traced back to the deviation of the 
velocity distribution from the Maxwell-Boltzmann one, 
w(v\x). At equilibrium, the velocity statistics is indeed 
given by 



p(v\x) = w(v\x) ( 1 



(5T(aO 



67T(x) 2 



-v • VT(aj) (15) 



except for terms of order 7~ 2 or higher. This deviation, 
albeit small, plays a crucial role as it breaks the velocity 
reversal symmetry along the temperature gradient direc- 
tion, and is therefore responsible for the anomalous en- 
tropy production. 

Further insight on the interpretation of the anoma- 
lous entropy can be gained by the following argument 
(see Figure [J. At the microscopic level (scales of or- 
der of T x / 2 /7, times of order I/7) the particle trajec- 
tory often "goes round in circles". Indeed, as seen from 



eq. (15), the trajectories with higher probability (posi- 



tive sign of the correction term) are characterized by a 
downgradient sweep with speeds larger than the thermal 
velocity, followed by a slow upgradient motion (red tra- 
jectory in Figure [J. Note that the correction does not 
contribute to the spatial flux of particles as can be di- 
rectly checked by multiplying by v and integrating over 
velocity. It does not contribute to the mean kinetic heat 
exchange either, since it is odd in velocity coordinates. 
During this futile cycling, however, heat is absorbed at 
higher temperatures and released at smaller ones, thereby 
producing entropy. In each cycle of duration ~ 7 _1 , 
the amount of entropy produced is approximately given 
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by the heat exchanged ~ T, times the difference of in- 
verse temperatures at the extremes of the trajectory, 
~ r- 2 |VT|(T 1 / 2 /7). Reversed cycles (such as the one 
depicted in blue in Figure [| give a similar contribution 
with opposite sign. However, forward cycles are more 
likely with probability ~ T -1 / 2 |VT|/7 as follows from 



eq. (15). This results in an overall average entropy pro- 
duction rate - 7- 1 T- 1 |VT| 2 that we recognize as the 
anomalous term in eq. (10). 



The present findings have a much broader range of ap- 
plicability than the simple example discussed here (see 
Supplementary material). The addition of potential and 
nonconservative forces is straightforward since the re- 
lated heat, work and entropy contributions all have reg- 
ular limits and are thus correctly described by the over- 
damped dynamics. A position-dependent friction, which 
may account for the temperature dependence of viscosity 
or hydrodynamic effects due to the presence of material 
boundaries, can also be included with minor changes. Fi- 
nally, it is also possible to consider time-dependent tem- 
perature, friction and forces, provided the variation is 
not faster than the timescale of overdamped motion. All 
these modifications do not alter our results. The inter- 
pretation is left unchanged as well. 

In conclusion, we have shown that thermodynamics at 
the micro-scale is intrinsically anomalous in the limit of 
vanishingly small inertia. The anomaly arises from the 
breaking of time-reversal symmetry that emerges when 
the particle is subject to a temperature gradient. How 
may the entropic anomaly impact observable phenom- 
ena ? The most conspicuous effect should be on the effi- 
ciency of thermal stochastic engines in nonuniform tem- 
perature environments. Indeed, the intrinsic irreversibil- 
ity arising from the anomalous contribution - which does 
not vanish even for quasistatic transformations - should 
irremediably hamper the ability of the engine of convert- 
ing absorbed heat into work. Recent advances in ex- 
perimental techniques of particle confinement, tracking, 
and heating might pave the way toward an experimental 
measurement of this effect [5H7]. 
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